function [ LL ] = loglikelihood(theta, id)
load('mainfile1.mat');
theta2=theta;
id2=id;
mis=mis(theta2, id2);
phi2=phi(theta2, id2);
mi=mi(theta2, id2, phi2);
pis=exp(mis)/(mi);
end

